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The most widely used approach for studying in the con- 
tinuum QCD effects that lie beyond the realm of pertur- 
bation theory are the Schwinger-Dyson (SD) equations. 
This infinite system of coupled non-linear integral equa- 
tions for all Green's functions of the theory is inherently 
non-perturbative and can accommodate phenomena such 
as chiral symmetry breaking and dynamical mass gener- 
ation. In practice one is of course severely limited in 
their use, and various approximations have been imple- 
mented throughout the years. Devising a self-consistent 
truncation scheme for the SD series is far from trivial. 
The main problem in this context is that the SD equa- 
tions are built out of unphysical Green's functions; thus, 
the extraction of reliable physical information depends 
crucially on delicate all-order cancellations, which may 
be inadvertently distorted in the process of the trun- 
cation. In order to partially compensate for this type 
of shortcomings, one usually attempts to supplement as 
much independent information as possible, by "solving" 
the complicated Slavnov- Taylor identities (STI), or by 
combining with results from lattice simulations. 

The truncation scheme based on the pinch technique 
(PT) 1, 2] implements a drastic modification already at 
the level of the building blocks of the SD series, namely 
the off-shell Green's functions themselves. The PT is 
a well-defined algorithm that exploits systematically the 
symmetries built into physical observables, such as S- 
matrix elements, in order to construct new, effective 
Green's functions endowed with very special properties. 
Most importantly, they are independent of the gauge- 
fixing parameter, and satisfy naive (ghost-free, QED-like) 
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Ward identities (WI) instead of the usual STI. The up- 
shot of this program is to first trade the conventional SD 
series for another, written in terms of these new Green's 
functions, and subsequently truncate it, keeping only a 
few terms in a "dressed-loop" expansion, maintaining at 
the same time exact gauge-invariance. 

Of central importance in this context is the connec- 
tion between the PT and the Background Field Method 
(BFM). The latter is a special gauge- fixing procedure 
that preserves the symmetry of the action under ordinary 
gauge transformations with respect to the background 
(classical) gauge field A®, while the quantum gauge fields 
A° appearing in the loops transform homogeneously un- 
der the gauge group As a result, the background 
n-point functions satisfy QED-like all-order Wis. The 
connection between PT and BFM, known to persist to all 
orders, affirms that the (gauge-independent) PT effective 
n-point functions coincide with the (gauge-dependent) 
BFM n-point functions provided that the latter are com- 
puted in the Feynman gauge ;4]. 

In this talk we consider the all-order diagrammatic 
structure of the effective gluon self-energy, H fiU (q), ob- 
tained within the PT-BFM framework. We explain that, 
as a consequence of all-order WI satisfied by the full 
vertices appearing in the corresponding diagrams, the 
transversality of ILj„ (q) is realized in a very special way: 
the contributions of gluonic and ghost loops are sepa- 
rately transverse. In particular, we study a truncated 
version of this new series, keeping only the terms of the 
gluonic one-loop dressed expansion, while still maintain 
exact gauge invariance. Our attention will focus on a de- 
tailed scrutiny of the necessary conditions for obtaining 
infrared finite solutions from the SD equation. 

Let us first define some basic quantities. First of all, 
it should be clear from the beginning that there are two 



1 



different gluon propagators appearing in this probiem, 
A l j, l ,(q), denoting the background gluon propagator, and 
A /JJ/ (g), denoting the quantum gluon propagator appear- 
ing inside the loops. In the Feynman gauge, A^fa) is 
given by 



P^(<?)A(q 2 



where the transversal projector, 
IV 0?) = a^v - 



% 1u 



(1) 



(2) 



The scalar function A(q 2 ) is related to the all-order gluon 
self-energy n^(q) by 



n M ,((7)=P^(g)n(g 2 ); A(q 2 ) 



iU(q 2 



(3) 



Exactly analogous definitions relate A /XI/ (q) with Tl llv (q). 

The diagrammatic representation of A~^{q) is shown 
in Fig.QJ. Notice that diagrams (62), (di), (cfe), (cfa) and 
(di), are characteristic to the BFM; within the PT they 
are generated dynamically, from the STIs triggered by 
the pinching momenta. 
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FIG. 1: The SD equation for the gluon propagator. Wavy 
lines with grey blobs represent full-quantum gluon propaga- 
tors, while the dashed lines with grey blobs denote full-ghost 
propagators. All external wavy lines (ending with a vertical 
line) are background gluons. The black dots are the tree-level 
vertices in the BFM, while black blob represents the full con- 
ventional vertices. The white blobs denote three or four-gluon 
vertices with one external background leg. 

As is widely known, in the conventional formalism the 
inclusion of ghosts is instrumental for the transversality 
of II™, (q) , already at the level of the one- loop calculation. 
On the other hand, in the PT-BFM formalism, due to 
new Feynman rules for the vertices, the one-loop gluon 
and ghost contribution are individually transverse |3(- 

As has been shown in ||, this crucial feature persists 
at the non-perturbative level, as a consequence of the 
simple Wis satisfied by the full vertices appearing in the 
diagrams of Fig.Q. Specifically, the gluonic and ghost 



sector are separately transverse, within each individual 
order in the dressed-loop expansion 

We will show this property for the one-loop dressed 
terms. We start by writing down the fundamental all- 
order WI for the full three-gluon vertex with one back- 

~ abc 

ground gluon, F MQ/3 , and for the full background gluon- 



ghost vertex F^ 

- abc 



- acb 



11^^/3(11,12^3) 
u — acb 

ii T tj. (92,91, 93) 



fjf 



abc 



Kl(i2)-Kl^) 



9f abC [D- 1 (q2)-D- 1 (q 3 )] ,(4) 



where on the RHS we have differences of inverse of the 
quantum gluon, A flL ,(q), and ghost, D(q), propagators. 

The closed expressions corresponding to the gluonic 
sector, at one-loop dressed expansion, (see Fig.Q) are 
given by 



9£(9)| 
n#(«)| 



1 f[dk] T^AZ(k)v b :' p x a 'a£,(* + q) 



5 [dk}T^A^(k) 



(5) 



with T aex a and T abex „ being the three and four bare gluon 
vertices in the Feynman gauge of the BFM [3j . 
For the ghost sector, we have 

H#(?)| bl = - f[dk}f^D ee ,(k)Tf X 'D xx ,(k + q), 



U ab (a)\ 



[dk] K b ;*D ex (k) 



(6) 



where T a ^ x and V ab ^ x represent the tree-level ghost-gluon 
vertices with one (two) background gluon(s) respectively 
0; the measure [dk] = d d k/(2ir) d with d = 4 - e the 
dimension of space-time. Observe that in our notation 
all the three and four-point functions with a tilde are 
vertices with at least one external (background) gluon 
leg. 

With the above WI we can prove that the groups (a) 
and (b) are independently transverse. We start with 
group (a) 



C A g'6 ab q^ [dk]M{k) 



l"K(l)L = -CA9 2 S ab q» j[dk]A"Jk), 



and thus 



(n£(ff)| 



(7) 



(8) 



Similarly, the one-loop-dressed ghost contributions 
give upon contraction 



2C A g l 8 ab q v [dk]D{k), 



9"n#(g)| = -2C A g 2 5 ab q„ j[dk]D(k), (9) 
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and so 



f (n»| bl 



nl(5)| 



o. 



(10) 



The proof of the individual transversahty of the groups 
(c) and (d), constituting the two- loop dressed expansion 
of H fJll/ (q), is slightly more cumbersome but essentially 
straightforward |5j. 

The importance of this transversahty property in the 
context of SD equation is that it allows for a meaning- 
ful first approximation: instead of the system of coupled 
equations involving gluon and ghost propagators, one 
may consider only the subset containing gluons, with- 
out compromising the crucial property of transversahty. 
More generally, one can envisage a systematic dressed 
loop expansion, maintaining transversahty manifest at 
every level of approximation. This is not to say, of course, 
that we have some a-priori guarantee that the subset of 
diagrams considered here is numerically dominant. Ac- 
tually, as has been argued in a series of SD studies, in the 
context of the conventional Landau gauge it is the ghost 
sector that furnishes in fact the leading contribution 6] 
Clearly, it is plausible that this characteristic feature may 
persist within the PT-BFM scheme as well, and we will 
explore this crucial issue in the near future 

Thus, in this formalism, the first non-trivial approx- 
imation for A _1 (<7 2 ) that preserves its transversahty is 
given by the gluonic terms of the one-loop expansion (di- 
agrams (di) and (02) in Fig.Q), written in closed form 
in Eq. ©. 

However, the equation given in (|SJ| is not a genuine 
SD equation, in the sense that it does not involve the 
unknown quantity A on both sides; instead, in the inte- 
grals of the RHS appears A. Replacing A by A is a highly 
non-trivial proposition, whose self-consistency is still an 
open issue. Its implementation may be systematized by 
resorting to a set of crucial identities relating A and A 
by means of a set of auxiliary Green's functions involving 
anti-fields and background sources. At this point we will 
assume that to a first approximation one may neglect the 
effects of the aforementioned auxiliary Green's functions, 
and carry out the substitution A — > A on the RHS of JSJ. 
Hence, the SD equation we will solve is written as 



A-V) 



q 2 +i 



& ai \q 2 )+& a ^ 



(11) 



where 



n£>(g) - \c A g 2 J[dk]r^A aa '(k)T, a/p/ A^'(k + q ) } 
-C A g 2 g^ f[dk]A(k). (12) 



(ii) introduce massless poles into the Ansatz for the full 
three-gluon vertex. 

The necessity of the first condition can be appreciated 
by observing that, on dimensional grounds, the value of 
A _1 (0) can only be proportional to two types of seagull- 
like contributions, 



n (a2) 



With the above equation at hand, we proceed to es- 
tablish the conditions necessary for obtaining infrared 
finite solutions for A _1 (q 2 ), i.e. solutions for which 
A _1 (0) ^ 0. There are two such conditions: one must 
(i) allow for non-vanishing seagull-like contribution, and 



% 



[dk] A(k), 



[dk]k 2 A 2 (k), (13) 



since, inside the diagrams of Ea. (|12fl . there can be at 
most two full gluon self-energies, A. However, it is well 
known that, due to the dimensional regularization rules, 
such contributions vanish perturbatively, ensuring the 
masslessness of the gluon order by order in perturba- 
tion theory. In order for finite solutions to emerge, one 
must assume that seagull-like contributions, such as those 
of Ea. l|13fl . do not vanish non-perturbatively. Naturally, 
this last step will force us to deal with the quadratic 
divergences, present in both integrals of Eci. lfLUl . and 
therefore a suitable regularization scheme must be sub- 
sequently employed. 

Allowing the non- vanishing of seagull-like terms is not 
the whole story however; one must determine in addi- 
tion the mechanism that will produce their appearance. 
One thing is certain: the seagull contributions determin- 
ing A _1 (0) do not originate from diagram (02) in Fig.Q. 
Instead, the required seagull contributions will stem from 
diagram (ai), after the inclusion of massless poles into the 
Ansatz for the full three-gluon vertex \T . Diagram (02) 
plays of course a crucial role in enforcing the transversal- 
ity of n M „ non-perturbatively, but in the absence of mass- 
less poles in the vertex one would still get A _1 (0) = 0. 

There is a relatively simple argument that amply 
demonstrates the subtle interplay between both require- 
ments. Specifically, iUfy (q) can be written in the gen- 
eral form 



m$\Q) = Q 2 Ml 2 )9^ + B(q 2 )q^ 



(14) 



where A(q 2 ) and B(q 2 ) are arbitrary dimensionless func- 
tions, whose precise expressions depend on the de- 
tails of the Fva'p' employed. The transversahty of 



(a 2 ) 
/if 



implies immediately the condition 



q 2 [A(q 2 ) + B(q 2 )] =lC A g 2 T 0l (15) 
and thus the sum of the two graphs reads 



-q 2 B(q 2 )P^(q). (16) 



Clearly from Ea. lfTH) . we conclude that A 1 (0) = 
lim (-q 2 B(q 2 )). 

Interestingly enough, once the transversahty of n AI „ 
has been enforced, the value of A _1 (0) is determined 
solely by lim (-q 2 B(q 2 )) . Evidently, if B(q 2 ) does not 

q 2 ^0 
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contain (l/<? 2 ) terms, one has that lim (—q 2 B(q 2 )) = 0, 

and therefore A _1 (0) — 0, despite the fact that To has 
been assumed to be non- vanishing. Thus, if the full three- 
gluon vertex T satisfies the WI of (@J), but does not con- 
tain poles, then the seagull contribution Tq ^ of graph 
(02) will cancel exactly against analogous contributions 
contained in graph (ai), forcing A -1 (0) = 0. 

We next proceed to study the SD of Eq. (fTTJ) . We will 
follow the methodology developed in [| and linearize 
the equation by resorting to the Lehmann representa- 
tion, together with a gauge-technique inspired Ansatz 
for the vertex J^ va p. This approximation yields a more 
tractable form for the resulting SD equation, which for 
the purposes of this preliminary analysis should suffice; 
of course, a non-linear study must eventually be carried 
out, and lies within our immediate plans. 

To simplify the form of the vertex required, we drop 
the longitudinal parts of A Mt/ inside the integrals, using 
A^ u (k) = —ig llu A(k). Omitting these terms does not 
interfere with the transversality of the external A^(q) 
5] . Then we obtain 



with 



\c A g 2 J [dk] Tf A(k)T uafi A(k + q) 



C A g 2 dg^ [dk]A(k), 



(17) 



= [2k + q) fl g a (3 - 2q a g^0 + 2qpg fia , (18) 



and 



q v T vaf} = A- l {k + q)-A- 1 {k) g af} . (19) 



The Lehmann representation for the scalar part of the 
gluon propagator reads 



A(9 2 ) 



dX 1 



p(X 2 ) 



q 2 — X 2 + ie ' 



(20) 



with no special assumptions on the form of the spectral 
density. 

This way of writing A(q 2 ) allows for a relatively simple 
gauge-technique Ansatz for T^ aj3 , which linearizes the 
resulting SDE. In particular, setting on the first integral 
of the RHS of Eq.|T7|) 



A(k) f va , A(k + q) = JdX 2 P (A 2 )^^- f U {k + q? \ x , + ie > (2D 
where T^ a p must be such as to satisfy the tree-level WI 

<ff^ - [(* + q) 2 - k 2 ] g aP = [(A + qf - A 2 ] g a p - (k 2 - X 2 )g a p . (22) 

Then, it is straightforward to show by contracting both sides of l|21|) with q", and employing (|20|) and (|22|l . that T^ Q/3 

satisfies the all-order WI of Eq. i|19[l . Of course, choosing T^ a g = r„ Qj a solves the WI, but as we will see in detail in 
what follows, due to the absence of pole terms, it does not allow for mass generation, in accordance with our previous 
discussion. Instead we propose the following form for the vertex 

PL/3 = f„ Q /3 + ci ((2k + q), + ^ [k 2 ~(k + q) 2 ] ) g a p + ^c 3 + [(k + q) 2 + k 2 ] ^ (q p g ua - q a9v p) , (23) 

which, due to the presence of the massless pole is expected to allow the possibility of infrared finite solution. Further- 
more, we treat the constants c\, ci and C3 as arbitrary parameters, in order to check quantitatively the sensitivity 
of the obtained solutions on the specific details of the form of the vertex. Notice that all new terms contributing to 
^va8 nave the correct properties under Bose symmetry. 

Thus, the vertex F entering in Eq. l|12[l .can be obtained as a combination of Eas. (l21|) and (|23|l . After rather lengthy 
algebraic manipulations of Ea. (|llfl (see |5j] for details), fixing C3 = ^ci to recover the perturbative result, we obtain 
for the renormalized A(g 2 ) (in the Euclidean space) 

A- 1 (q 2 ) = q 2 { K + bg 2 \ dz[l- — ) A(z) \ + jbg 2 \ dzz[l--^) A(z) + A" 1 (0) , (24) 



r 



with b = 10CA/487T 2 , and 

ba 2 <T r ~ - 6(ci+c 2 ) _ 4 + 4ci+3c 2 
A-\0) = - b -^Jd*kA(k 2 ), (25) CT = 5 ' 7= 5 ■ (26) 
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The renormalization constant K is to be fixed by the 
condition, A _1 (/_t 2 ) = fi 2 , with fj, 2 > A 2 . Notice that 
the deviation of b from the value b = 1ICU/487T 2 , the 
standard coefficient of the one-loop (3 function of QCD, 
is due to the omission of the ghosts. From Ij25(l it is clear 
that when a = automatically A _1 (0) vanishes, despite 
the inclusion of (02). Note however that having poles is 
not a sufficient condition: if c\ — — C2, there is no effect. 

It is interesting to study the UV behavior for A(q 2 ) 
predicted by the integral equation (|24p. At large q 2 we 
can safely replace the factors (I — 4z/q 2 ) 1 / 2 — > 1, arriving 
at the following simplified version of that equation, 



l + bg 2 j dzA{z)\ , (27) 



whose solution can be easily obtained by casting it into an 
differential equation, written in terms of the form factor 
G{q 2 ) = q 2 A(q 2 ), which lead us to 



A-\q 2 ) = q 2 



1 + 2bg 2 In 



1/2 



(28) 



Obviously, upon expansion this expression recovers the 
one-loop result A~ 1 (q 2 )\ pert = q 2 ^1 + bg 2 \n(q 2 / y 2 )^ 

correctly, but A(q 2 ) does not display the expected RG 



behavior at higher order. The fundamental reason for 
this discrepancy can be essentially traced back to having 
carried out the renormalization subtractively instead of 
multiplicatively 0, , a fact that distorts the RG struc- 
ture of the equation. 

As is well-known, due to the Abelian WI satisfied by 
the PT effective Green's functions, A^ 1 (q 2 ) absorbs all 
the RG-logs, exactly as happens in QED with the photon 
self-energy. Consequently, the product d(q 2 ) = g 2 A(q 2 ) 
should form a RG-invariant (^-independent) quantity. 
Notice however that Ea. (|24fl does not encode the cor- 
rect RG behavior: when written in terms of the RG in- 
variant quantity d(q 2 ) = g 2 A(q 2 ) it is not manifestly 
<7 2 -independent, as it should. 

In order to restore the correct RG behavior at the level 
of (|24(l . observe that such equation requires an extra 
power of g 2 in their integrands on the RHS. Then, we 
use the simple prescription whereby we substitute every 
A(z) appearing on RHS of Ea.pi) by 0,0 



A(z) 



g 2 A(z) 



9 2 {z) 



= [l + bg 2 Hz/iJL 2 ))A(z). (29) 



which allows us to cast Eq.J23J) in terms of the RG- 
invariant quantities d{q 2 ) and g 2 (q 2 ) in the following way: 



d-\q 2 ) = 




dz 1 - 



— / dz z 



4z 


V /2 




/ 9 2 (z) / 


Az \ 


1/2 d(z) \ 


V 2 ) 


9 2 (z) 



d-\o) , 



(30) 



where 



(31) 



It is easy to see now that Ea. (l3()ll yields the correct UV 
behavior, i.e. d~ x {q 2 ) = b q 2 \n(q 2 / A 2 ) . 

When solving (|3f)|> we will be interested in solutions 
that are qualitatively of the general form 



d(q 2 ) = 



9lAl 2 ) 



m 2 (q 2 



where 



—2 



b In 



g 2 + /( g 2 ,m 2 (g 2 )) 
A 2 



(32) 



(33) 



The quantity g 2 ^^ 2 ) represents a non-perturbative ver- 
sion of the RG-invariant effective charge of QCD: in the 



deep UV it goes over to g~ 2 (z), while in the deep IR it will 
be finite due to the presence of the function f(q 2 , m 2 (q 2 )), 
whose form will be determined by fitting the numerical 
solution. 

The function m 2 (q 2 ) may be interpreted as a momen- 
tum dependent "mass" . On general arguments dynam- 
ically generated masses must vanish asymptotically. In 
order to determine the asymptotic behavior that Eq. (f 3()|) 
predicts for m 2 (q 2 ) at large q 2 , we replace Ea. H32(l on 
both sides of Eq.!^, set (I - Az/q 2 ) 1 ^ 2 -> 1, and de- 
mand the consistency of both sides, obtaining finally 



m\q 2 ) 



„,2 1 -a 

m n m 



g 2 /A 2 



with 



1 + 7 • (34) 



Indeed the gluon mass vanishes at UV as an inverse power 
of ln(g 2 ), since a > 0. Actually, as we will see below, 
the regularization of Ea. pil) imposes a more stringent 
constraint, requiring that 7 > 0, thus restricting through 
Ea. (|26|) the possible values of c\ and C2 in ^ a g- 
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As mentioned before, the seagull-like contribution (de- 
noted collectively by d~ 1 (Q) in Eq.(j20J) are essential for 
obtaining IR finite solution for d(q 2 ). However, the in- 
tegral (|31|) should be properly regularized, in order to 
ensure the finiteness of such a mass term. 

For the regulation of the quadratic divergences present 
in the integral (|31|l . we rely on two basic ingredients: (i) 
the standard integration rules of the dimensional regu- 
larization and (ii) a constraint on the allowed values of 
the anomalous mass-dimension a. 

With this in mind, we recall that according to the di- 
mensional regularization rules, J[dk]/k 2 = 0, allowing us 
to rewrite the Ea. (f3"TJ) (using ffi^i ) as 



ba 



[dk] 



1 

[k 2 + m 2 {k 2 )]g 2 (k 2 ) ~ k 2 
m 2 (k 2 ) 



k 2 [k 2 + m 2 (k 2 )] 

f(k 2 ,m 2 (k 2 )) 



~ /[(//,•] .71 /,--! In ( I 



(35) 



The inspection of the two integrals on the RHS sep- 
arately reveals that, if m 2 (k 2 ) falls asymptotically as 
power of h\~ a (k 2 ), with a > 1, then the first integral 
would converge, by virtue of the elementary result 



dz 



1 



z (In z 



|l+7 



7 (In z)"' 



(36) 



which, of course, requires that 7 > 0. The second integral 
will converge as well, provided that f(k 2 ,m 2 (k 2 )) drops 
asymptotically at least as fast as ln~ c (fc 2 ), with c > 0. 
If, for example, / = pm 2 (k 2 ) (with 1 + 7 > 1, for the first 
integral to converge), then the convergence condition for 
the second integral is automatically fulfilled. Notice that 
perturbatively <i~ 1 (0) vanishes; this is because m 2 (k 2 ) = 
to all orders, and therefore, since in that case also / = 0, 
both integrals on the RHS of l|35[) vanish. 

Evidently, Eas. H30(l and i|35|l form a system of equa- 
tion; the role of the first is to provide a solution for the 
unknown RGI quantity, d(q 2 ), while the second acts as 
an additional constraint, restricting the number of pos- 
sible solutions. Therefore, Eas. (|3T7|) and (|33|l should be 
solved simultaneously. 

Using an iterative method, we performed a detailed 
study of these two equations, where for each d _1 (0) cho- 
sen we vary 7 and a in order to scan the two-parameter 
space of solutions. 

In Fig.J2J, w e show numerical results for d(q 2 ), for dif- 
ferent values of cZ — 1 (0). All these solutions satisfy the 
constraint imposed by Ea. (|3"5j) and their respective val- 
ues for a are described in the inserted legend. As ex- 
pected, the gluon propagator behaves at low momenta 
like a constant, whose value is determined by <i(0); in ad- 
dition, it obeys the correct ultraviolet behavior. As can 
be observed from the plot, d(q 2 ) starts off as constant 
until the neighborhood of q 2 = 0.01 GeV 2 , where the 



curve bends downward in order to match with the per- 
turbative asymptotic behavior at a scale of a few GeV 2 . 
All these solutions can be perfectly fitted by Eq.(J22J); the 
functional form of f(q 2 ,m 2 (q 2 )) and m 2 (q 2 ) that better 
describe our data sets are given by 



f(q 2 ,m 2 (q 2 ))=p 1 m 2 (q 2 ) 
and the dynamical mass, 
m 2 (q 2 ) = '/ 



P2 



m A (q 2 ) 
q 2 + m 2 (q 2 ) 



In 



q 2 + p 1 m 2 



A 2 



In 



Pi ml 



A 2 



(37) 



, (38) 



with exponent a = 1 + 7. 
In all cases, we fixed pi 



= 4; therefore, the unique 



free parameter is P2, since the value of Wq, which is also 
related to p2, can be directly obtained by setting q 2 = 
in Ea . 1)32(1 . From this follows immediately that the value 
of the infrared fixed point of the running coupling, r(0) 
will be determined by the value assumed by /(0, ttiq), 



g- p 2 (0) = 61n 



(4 + p 2 )m 2 
A 2 



(39) 



Obviously, the maximum value obtained for g 2 p (0), is 
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FIG. 2: Results for d(q 2 ) fixing different values for d _1 (0) (all 
in GeV 2 ). All these solutions satisfy the condition given by 
Eq. 1351 . Their respective values for a are given in the legend, 
in all cases we set C2 = in Eq, l26H . 

the one that minimizes (4 + p 2 ) an d, at same time, 
keeps it bigger than A 2 , in order to avoid the pole. Of 
course, if we set p 2 — in Eq.(l37|l. we automatically re- 
cover the solution proposed in [l| ; however our numerical 
solution requires bigger values for the coupling ~g 2 ]P {q 2 ), 
and therefore pi assumes negative values, as can be ob- 
served on Fig.J2J) 

The running charges, a{q 2 ), for each solution presented 
on Fig.Q are displayed in Fig.©. Observe that as the 
value of <i -1 (0) decreases, the value of the infrared fixed 
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FIG. 3: The running charge, ct(q 2 ), corresponding to the 
gluon propagator of Fig.©. Clearly, a(0) increases as d _1 (0) 
decreases 
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FIG. 4: The ratio mo/A as function of the parameter a. The 
increase is exponential, given by mo/A = Ai exp(er/ti) + yo, 
where Ai = 0.775, h = 0.25 and y = 0.436. 

point of the running coupling, a(0), increases. Accord- 
ingly, from Fig.© and Fig.©, we can conclude that, if 
smaller values of a were to be favored by QCD, the freez- 
ing of the running coupling would occur at higher values. 
It should also be noted that the values of a(0) found here 
tend to be slightly more elevated compared to those of 
Q (for the same value of ttiq) . 



Finally, we analyze the dependence of the ratio mo/ A 
on a; the former is extracted from Eci. (|32|) by setting 
q 2 = 0. This dependence is shown in the Fig.©, corre- 
sponding to the cases presented in Fig.©. We observe 
that as we increase the value of cr, namely the sum of the 
coefficients of the massless pole terms appearing in the 
three gluon vertex, the ratio mo/ A grows exponentially. 

In conclusion, we have presented an analysis of the var- 
ious intertwined issues involved in the study of dynam- 
ical gluon mass generation through SD equations. The 
analysis was carried out in the context of the PT-BFM 
scheme, where various crucial properties are preserved 
manifestly. Most notably, the transversality of the non- 
perturbative gluon self-energy is enforced order by order 
in the dressed loop expansion and separately for gluons 
and ghost. We have seen in detail that for the existence 
of infrared finite solutions two requirements are indis- 
pensable: (i) the non-vanishing of seagull-like contribu- 
tions beyond perturbation theory, and (ii) the presence 
of massless poles in the trilinear gluon vertex. The result- 
ing equation was linearized by resorting to the Lehmann 
representation and a gauge-technique inspired solution 
of the corresponding WI. A simple Ansatz for the three- 
gluon vertex was constructed, which contains massless 
poles, thus allowing for the appearance of infrared finite 
solutions. This vertex satisfies the correct WI, but oth- 
erwise is purely phenomenological, in the sense that it is 
not QCD-derived, nor does it contain the most general 
Lorentz structure. Numerical solutions were obtained for 
the RG-invariant quantity d(q 2 ); they can be fitted per- 
fectly by means of a running coupling that freezes in the 
IR, and a dynamical mass that vanishes in the UV. We 
have found that the actual values of a(0) depend strongly 
on the combined strength of the pole terms appearing in 
the vertex. This strongly suggests that the value of this 
IR fixed point should be determined by means of a de- 
tailed non-perturbative study of the three-gluon vertex, 
either on the lattice or through its own SD equation. 
Needless to say, many of the issues considered in this 
talk are far from settled, and a lot of independent work 
is necessary before reaching definite conclusions. 
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